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Abstract 

We experimentally and theoretically study electromagnetic properties of optically thin silicon 
carbide (SiC) membranes perforated by an array of sub-wavelength holes. Giant absorption and 
transmission is found using Fourier Transformed Infrared (FTIR) microscopy and explained by 
introducing a frequency-dependent effective permittivity e^siuj) of the perforated film. The value 
of eeff(w) is determined by the excitation of two distinct types of hole resonances: a delocalized 
slow surface polariton (SSP) whose frequency is largely determined by the array period, and a 
localized surface polariton (LSP) which corresponds to the resonances of an isolated hole. Only 
SSPs are shown to modify eefr(w) strongly enough to cause giant transmission and absorption. 
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Diffraction of light is the major obstacle to increasing the density of optical circuits and 
integrating them with electronics Q|: light cannot be confined to dimensions much smaller 
than half of its wavelength A/2. Utilizing materials with a negative dielectric permittivity 
circumvents diffraction limit because interfaces between polaritonic (e < 0) and dielectric 
(e > 0) materials support quasi-electrostatic waves (surface polaritons) that can be confined 
to sub-A dimensions. Negative e can be due to either collective oscillation of conduction 
electron in metals (plasmons) 0| or lattice vibrations in polar crystals such as ZnSe, SiC, 
InP P|. Both metallic and SiC interfaces have been shown to guide surface plasmon (and 
phonon) polaritons over long distances 0.0. 01, support negative-index waves , and 

serve as elements of a super-lens jiol 111! Therefore, the explosion of research aims at 
understanding basic polaritonic components: grooves in metal films |13i, nanoparticle-film 



polaritons 



14| . hole arrays |1.5], and single holes [iG^ . 
Sub-A hole arrays were the first to attract significant attention after the discovery of 
extraordinary optical tra^n^i^ion M through perforated opt.cally thick rrretallic filn. and 

the ensuing controversy [121, Hq] about the role of surface polaritons in transmission enhance- 
ment. The difficulty with interpreting the role of surface polaritons (SPs) in transmission 
through optically thick (Askin = '^/4vrA/|e| ^ H , where H is the film thickness) films arises 
from the fact that even holes in a perfect conductor (which does not support surface po- 
laritons) can spoof surface polaritons especially when A is close to array period L. No 
such spoofing occurs for optically thin films. 

In this Letter we report results of the spectroscopic study of such suspended perforated 
SiC membranes, with both hole diameter D <^ X and period L < X significantly sub-A. SiC is 
chosen because of its low losses \^ and the possibility of growing high quality thin SiC films on 
Si substrates [2^. While the majority of surface phonon polariton experiments to-date have 
been conducted with SiC substrates 0,0, 0|; ^^^^ membranes hold a greater promise 



for fabricating novel multi-layer nanophononic structures [1^. Using FTIR microscopy, we 
demonstrate that, in adjacent frequency ranges, giant transmission and absorption of the 
incident radiation is realized. Both phenomena are shown to be due to the excitation of quasi- 
electrostatic SPs. It is shown that a perforated film can be described as a metamaterial with 
the effective permittivity e^s = + iei strongly modified by the excitation of SPs: regions 
of giant transmission and absorption are related, respectively, to the lowering of \er\ and 
increase of e^. In addition, we show theoretically that two types of SPs are supported by the 
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hole arrays in a polaritonic membrane: localized surface polaritons (LSPs) and delocalized 
SSPs. 

Square arrays of round holes were milled by a dual focused ion beam (FIB) system (FEI 
Strata 235) in a suspended SiC membrane. The starting material for membrane fabrication 



was a 458 nm thick single-crystalline 3C-SiC fi^ 



ished) side of a 0.5 mm thick Si(lOO) wafer [20!]. The 0.35 x 0.35 mm suspended membrane 



m heteroepitaxially grown on the front (pol- 



121. All 



was produced by anisotropic KOH etching of the Si from the back (unpolished) side 
holes were round in shape (see insets to Fig. [TJ, with the edge radius of curvature of order 
30 nm. The hole diameter and period in different samples varied in the range 1 < D < 2/im 
and 5 < L < 7/im, and a typical perforated area was 150 x 150/im in size. Inset of the 
Fig. n shows one fabricated sample with D = 2fim and L = 7fim. An FTIR microscope 
(Perkins-Elmer Spectrum GX Autolmage) with spatial domain (100/im)^ has been used to 
measure the transmission and reflection from the perforated and non-perforated areas in the 
1500 — 700cm~^ spectral range. The FTIR microscope utilized Cassegrain optics with the 
angle of incidence between 9° — 35°. Input IR radiation was polarized by a wire- mesh (ISP 
Optics) polarizer. To quantify the effect of the holes, transmission and absorption spectra of 
a non-perforated SiC membrane are subtracted from the corresponding hole array spectra 
and plotted in Fig. ^for several values of D and L. 

Below we focus on the {D = 2fim, L = 7fim) sample. Although the deeply sub-A holes 
occupy only rj = 6% of the total sample area, they increase the transmittance through 
the sample by Th = 20% (from Tgim = 10% to Tpcrf = 30%) at \tr = 11.95/xm. This 
x3 transmission increase over the fractional holes' area is quite remarkable; outside of the 
reststrahlen {e{uj) < 0) region of SiC the increase is much smaller than 77. Even more 
remarkable is the increase of absorption (up by almost 40% from < 1% in a non-perforated 
film) at Xabs = 11.8/im. Other perforated samples {L = 5/im, D = 1/im) and {L = 7/im, 
D = lyum) were also experimentally studied and gave qualitatively similar results as shown 
in Fig.^ The absorption maximum is blue-shifted for smaller holes {Xabs = 11.57/im for the 
L = 7^m, D = 1/im sample) and for smaller periods {Xabs = 11.49/im for the L = 5/im, D = 
1/im sample), in accordance with the theory presented below which explains the anomalies 
in transmission and absorption as manifestations of surface phonon polariton excitation. 
This theory paves the way to the engineering of optical properties of polaritonic films using 
resonant excitation of surface polaritons. More complex three-dimensional materials can 
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then be built using multi-layer polaritonic films as building blocks. 

Because of the sub-A nature of the perforated film, it can be described as an effective 
medium with a frequency-dependent permittivity eefr(i^) determined by several strongest 
electrostatic (ES) resonances. Strengths and frequencies of these resonances are numeri- 
cally computed u sing two different formalisms: generalized eigenvalue differential equation 
(CEDE) 
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2J| and the surface integral eigenvalue equation 



26|. The former 



method has been successfully applied to periodic sub-A polaritonic crystals consisting of 
non-connected polaritonic inclusions 2^. To the best of our knowledge, this is the first 
generalization of this method to the system with a continuous polaritonic phase. 

The details of the GEDE approach as applied to periodic nanostructures consisting of two 
material components (one with a frequency-dependent e{uj) < and another with = 1) 
is described elsewhere js^El- Briefly, the following steps are followed. First, the GEDE 



V ■ 



= SjV^0j is solved for the real eigenvalue s,, where 6{x) = 1 inside the 
polaritonic material and 6{x) = elsewhere, and (pi is a potential eigenfunction with the 
dipole symmetry periodic with period L in the plane of the film. Second, the dipole strengths 
fi proportional to the squared dipole moment of each resonance are calculated [2^ • Finally, 
the ES permittivity eefr(co') is found by summing up the contributions of all dipole-active 
resonances: 
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f 



siuj) ^ siuj) — Si 



(1) 



where s = [1 — e{uj)Y^ serves as a frequency label. Note the pole at s = (corresponding to 
e —00) in Eq. (Q) which emerges due to continuous polaritonic phase. "Frequencies" Sj 
and strengths /, of the dominant dipole resonances computed for the experimentally relevant 
parameters {D = 2fim, H = 458nm, and L = 7/im) are: Sq = (eo = — c>o) with /o ~ 0.88, 
si ^ 0.1241 (ei ^ -7.058) with /i ^ 0.041, S2 ^ 0.1958 (es ^ -4.107) with /a ^ 0.0054, 
S3 ^ 0.255 (63 ^ -2.922) with /s ^ 0.0036, and S4 ^ 0.666 (64 ^ -0.50) with f^ ^ 0.0049. 
Other hole diameters were also studied, and the corresponding €1^2 are indicated by squares 
in the inset of Fig. 01 Note that eefr(co') is always finite because s is a complex number due 
to finite losses in SiC while all Si are real. 

Experimental results depicted in Fig. ^ can now be interpreted using ees{^)'- the normal 
incidence transmission and reflection coefficients through a slab of thickness H are given by 
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where n = y/e^,kn = u/c and ri = (1 — + n). The standard polaritonic formula 

recently verified for suspended SiC membranes was used for esic- We have calculated 



transmittance T and absorbance A = 1 — R — T of the perforated film, subtracted the 
corresponding quantities for the non-perforated film (eefr replaced by esic), and plotted the 
corresponding differential quantities in Fig. |21 To facilitate the interpretation of the giant 
transmission and absorption effects, a segment of Re[eefr], Ini[eefr] dependence is plotted in 
the inset to Fig. linear the strongest resonance at Ai = 11.3/im, where Re[esic(Ai)] = ei. 

The absorption spike (dotted line) is due to the peak of Im[ecfr] at A = Ai. The trans- 
mission maximum occurs due to the decrease of the absolute value of Re [eefr] at A = Xmax 
shown in the inset (solid line) to Fig. |21 Note that the enhanced transmission predicted from 
eefj occurs for Xmax > Ai in agreement with experimental observations. Other much smaller 
absorption and transmission peaks are also predicted. However, experimental accuracy was 
not sufficient to identify these additional variations of T and A in Fig. ^ Because spacing 
between the holes {L = 7/im) was not negligible compared with the wavelength (A ~ 11/im), 
we have carried fully electromagnetic (EM) calculations of T and A using the finite-elements 
frequency domain (FEFD) solver COMSOL. The results displayed in Fig. El show qualita- 
tive agreement with the eog-based calculation, with the exception that all transmission and 
absorption maxima are slightly red-shifted. This red shift is a previously noted [2^ phe- 
nomenon explained by the EM corrections to the purely ES response of sub-A polaritonic 
structures. 

One manifestation of resonances of ecsi^^) is an increase of the peak electric field -Emax 
inside the hole in response to the applied across the structure ac electric field with amplitude 
Eq. The ratio of E^g_^/ Eq is plotted in Fig. [Has a function of A for a fixed hole diameter 
D = 2/im and variable period L. Four enhancement spikes corresponding to resonances of 
eefr(co') can be identified for all periods. Frequencies of the three red resonances (A > 10.5/im) 
are located in the Re[esic] < band and are all strongly dependent on L. On the contrary, 
blue resonance (Aioc ~ 10.45/im) belongs to the —1 < Re[esic] < range and is period- 
independent. Another striking difference between the red and blue resonances is that the 
former are very de-localized, while the latter is strongly localized near the hole (see the two 
insets to Fig. Hj). Identification of the period-independent resonance as the LSP of a single 



hole, and of the period-dependent resonances as delocahzed waves related to SSPs of the 
smooth film is one of the main theoretical advances of this Letter. 

To validate this identification, we start with reviewing the properties of SPs of a smooth 
thin negative-e film surrounded by vacuum. Through symmetry, SPs can be labeled by 
the parity of the in-plane electric field with respect to the mid-plane, and by a continuous 
in-plane wavenumber k. A dispersion relation u vs. k in the ES limit of k ^ u/c is given 
implicitly by €±{00) = — tanh(/ci?/2)^^, where ± refers to even and odd modes, respectively. 
Note that even ("slow") modes are located in the —00 < e < — 1 part of the spectrum, 
while the odd ("fast") modes are in the — 1 < e < range. EM corrections to the dispersion 
relation do not alter this ordering (except ioi k ^ u/c case that is not relevant to this work). 
The EM dispersion curve and its ES approximation (expressed as Re[esic](co') vs. k) for the 
fast and slow SPs on a smooth SiC film are plotted in Fig. 01 

Because k > uj/c for SPs, they cannot be excited by an EM wave incident from vacuum. A 
perforated film, however, acts as a crossed diffraction grating and provides coupling between 
normally incident radiation (with almost uniform electric field Eq = EqCx) and surface modes 
of the film with the in-plane wavenumbers k(^rn,n) = \'^{jnx + riy)\ = "^^/m? + j^]. 
Resonances of egg (tu) can be understood as the result of strong coupling of Eq to slow surface 
polaritons (SSP(m,n) ^28|). This interpretation is verified by extrapolating the diameter- 
dependent resonant e for the two low-frequency eigenmodes to D = shown in the inset 
to Fig. El It is apparent from Fig. (HI that these two resonances correspond to excitation of 
SSP(0, 1) and SSP(1, 1) of the smooth film. Another conclusion that can be drawn from 
Fig. El is that EM effects shift the frequencies of SPs towards red. This explains the red shift 
of transmission and absorption spikes in EM simulations from their positions in ES model. 

No identification with one of the smooth film SSPs can be made for the LSP resonance at 
A = 10.45/im. In fact, this highly localized near the hole perimeter surface wave resonance 
can be thought of as an even-parity "defect state" created by the presence of a single hole 
in a negative-e film. Because the frequency range for which —1 < e{uj) < is a stop- 
band for even-parity propagating SPs, the even-parity LSP can exist in this, and only this 
frequency range. Because of the localized nature of the LSP, its frequency is insensitive to 
the proximity of other holes (i. e. to the period L) but is sensitive to the aspect ratio of the 
hole and even to the radius of curvature of hole edge. 

In conclusion, we have experimentally demonstrated giant absorption and transmission 
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through optically thin SiC membranes perforated by a rectangular array of round holes. 
These effects are theoretically explained by introducing the effective permittivity of a perfo- 
rated membrane with a complex frequency dependence due to resonant coupling to almost- 
ES surface phonon polaritons. Two types of phonon polariton resonances are theoretically 
uncovered: (i) delocalized modes related to the SSPs of a smooth film, and (ii) a localized 
surface polariton (LSP) of a single hole that exists in the spectral range complementary to 
that of SSPs. 

This work is supported by the ARO MURI W911NF-04-01-0203 and the DARPA contract 
HR0011-05-C-0068. We gratefully acknowledge Drs. A. Ahev and A. A. Zakhidov for their 
assistance with FTIR microscopy and Dr. C. Zorman for growing SiC films and useful 
discussions. 
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FIG. 1: (Color online) Extra transmission (blue) and absorption (red) of square arrays of circular 
holes in a 458nm thick SiC film, relative to non-perforated film. Solid: L = 7;um, D = 2^m] 
dash-dotted: L = Turn, D = 1/im; dotted: L = 5/xm, D = 1/xm. Left inset: SEM image of the 
L = 7fim, D = 2;um sample. Right inset: SEM image of a D = 2fj,m hole (at 52° from normal). 



9 





1 


1 


r5 1 1 




0.4 


p — p 1 
eff SiC 1 




M 




1 
1 




3 : ; ' 


1 
1 


0.3 


I 
1 

-vS/' 


1 
I 
1 

\ 


2 ; ; ' 

1 ; ' 










■ 1 




0.2 






-1 ; W 


AX 


0.1 



_ 10.5 11 

-s/|'V 


11.5 12 ' I - ' ' 




-0.1 






\J \ "\ 
\ 1 \ 
1 1 \ 1 




-0.2 






II \ / 
* ' V 





10.4 10.6 10.8 11 11.2 11.4 11.6 11.8 12 
wavelength, micron 



FIG. 2: (Color online) Theoretical extra transmission (green solid and dash-dotted curves) and 
absorption (red dashed and dotted curves) of a array of 2/xm holes in 458 nm film of SiC. Solid 
and dashed curves: FEFD simulation of EM wave scattering; dash-dotted and dotted: theoretic 
estimate based on ES egg, which is plotted on the inset. 




FIG. 3: (Color online) Dispersion relation of the even (green, lower curves) and odd (blue, upper 
curves) surface polaritons on a SiC film with thickness H = 458 nm. Solid lines: exact dispersion 
relation in the form of Re[esic('^)] vs. k, dashed lines: ES approximation. Two vertical lines: 
k = 2tt/L and k = 2-v/27r/L for L = Ijim.. Inset: position of delocalized resonances SSP(1,0) and 
SSP(1, 1) as a function of hole diameter D (squares on dashed line - ES eigenvalue simulations; 
stars on solid line - absorption peak positions in EM FEFD simulations). 
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FIG. 4: (Color online) Electric field enhancement in the symmetry plane of a SiC film perforated 
with a L = 7/im-period square array of D = 2/im round holes. Insets: ES potential profile at the 
resonances: (left) LSP resonance, and (right) SSP(1,0) resonance. 



11 



